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We analyze the mechanism that determines the boundary 
of stability in Taylor-Couette flow. By simple physical argu- 
ment we derive an analytic expression to approximate the 
stability line for all radius ratios and all speed ratios, for 
CO- and counterrotating cyhnders. The expression includes 
viscosity and so generalizes Rayleigh's criterion. We achieve 
agreement with linear stability theory and with experiments 
in the whole parameter space. Explicit formulae are given for 
limiting cases. 

PACS 47.15.Fe, 47.20.Gv 



1. Introduction 

The Taylor-Couette system has been studied for more 
than 100 years, beginning with the first experiments by 
Couette y , Mallock [§, and the pioneering work by Tay- 
lor ; for a recent review see Tagg [|| . 

We here focus on the question when and why the la- 
minar flow becomes unstable. This can and has been 
evaluated with linear stability theory I^ EI-fl , and was 
measured for a variety of radius ratios P]lC|-[l^ . The 
aim of this note is to elucidate the physical mechanism 
and derive a simple approximate analytical formula for 
the stability boundary in the whole parameter space. 

Earlier attempts of this kind have been made by Don- 
nelly and Fultz Coles @, Joseph §, and Eckhardt 
and Yao [|l5|. However, these authors did not consider 
all relevant aspects. Those include the smooth crossover 
from a solid boundary (for corotation) to a free one (for 
sufhciently strong counterrotation) and the localization 
of the initial perturbation in the center of the relevant 
gap. We elaborate that in Sects. 2-4. 

We consider an incompressible fluid with kinematic vis- 
cosity v flowing between two concentric cylinders with 
inner radius ri and outer radius r2, and angular veloci- 
ties LOi > and W2, positive for corotation and negative 
for counterrotation. This fluid flow is characterized by 
the Reynolds numbers Ri = uJiTid/v [i = 1,2), where 
d — r2 — ri is the gap si ze between the cylinders. 

The well-known 16| laminar flow is purely azimuthal 
and has the velocity profile 



Uy,{r) = Ar + 



B 



with 
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The radius ratio r/ = ri /r2 characterizes the geometry of 
the system. For counterrotation a nodal surface exists 



where 



0, 



r2 



I Rl-T]R2 
Rl-'q-^R2 



i?2 < 0. 
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2. Derivation of the stability criterion 

Let us start with Rayleigh's stability argument for in- 
viscid flow : Consider two neighbouring fluid rings at 
r and r + (5r of the laminar flow. Stability is guaranteed 
if their hypothetical exchange by virtual motion through 
the surrounding flow costs kinetic energy. Neglecting 
momentum loss by viscosity, the angular momenta per 



mass L{r) 



,{r) of the rings are conserved, and with 



their kinetic energy per mass being m^/2 we arrive at 
Rayleigh's famous stability criterion L^(r)' > when 
considering the energy balance 



kin 



L{r + Sr) 
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L{r + Sr) 
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L^ir)' 



{Srf > 0. 



5r has been choosen infinitesimal and the prime denotes 
the r-derivative. Stability is predicted for corotation if 
Ri < R2/TIJ ^-iid in the case of counterrotation for Ri — 
only. 

Now we take into account viscosity. During the motion 
to exchange the two fluid rings as described above they 
loose angular momentum and thus kinetic energy by dis- 
sipation. The frictional energy loss per mass is given by 
the viscous force per mass times the shift Sr, 

SE^is ~ (-^ Sur)Sr. 

£ denotes the typical diameter of the ring (in radial as well 
as in axial direction). Sur — Sr/6t is the radial velocity of 
the shift which takes place within time St. The sym- 
bol means 'up to a factor of order unity'. For instability, 
(5i?kin has to exceed (5i?vis- Its lower limit is related to the 
upper limit for St, which is set by the diffusive time scale, 
on which the rings loose angular momentum, St ~ i'^/v. 
We now obtain from (5i?kin ~ SE^is 
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as a relation which characterizes the stability boundary, 
or, after inserting the profil (0), i.e., L — Ar^ + i3. 
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Here = ~B/A, which for negative R2 coincides with 
the squared position of the nodal surface. 

In contrast to the inviscid case the timescale for the 
motion and the size of the rings enter here; the expression 
on the rhs. of has also been used in Q. As we will 
see soon, the radial position of the rings is important as 
well. 

£ is some small but finite fraction of the part of the 
gap which is susceptible to instability. For positive R2 
this is the whole gap between the cylinders, therefore we 
write £ = ad with a a constant considerably smaller than 
1 . We shall determine it by comparison with linear 
stability analysis. 

In the counterrotating case L^{r) decreases in the 
range from the inner cylinder up to the nodal surface 
only. We therefore compare £ with d„ = r„ — ri, which 
is an old idea, cf also [1^] . We further take into account 
that the flow patterns appearing at the instability line 
extend beyond the nodal surface . This may be under- 
stood as a change in the type of boundary from solid (for 
R2 > 0) to sort of free (for R2 < 0). The simplest way 
of describing this crossover is 



£^adA(a'^) with A(x) 
d 



if X <1 
if a; > 1 



(5) 



Thus, £ = aad„ for R2 < R2, where R2 is determined by 
the condition adn/d =1. If a > 1 we get R2 < 0. This 
delay in the decrease of £ when R2 becomes negative can 
also be seen in Fig. 12 of Ref. for the axial spacing of 
vortices. We expect this crossover to be even a smooth 
process. Therefore we demand that the stability curve 
Ri{R2) is at least continuously differentiable. This will 
fix the parameters a and R2 ■ Coles uses the same ansatz 
£ = aadn but for the whole range i?2 < 0, which 
makes the stability boundary discontinuous at i?2 = 0. 

Now we address the remaining question at which po- 
sition Tp the virtual motion of the rings has to be con- 
sidered. At first sight, one might think to choose it as to 
maximize the virtual energy gain. This leads to rp = ri. 
But more reasonable seems a position where the fluid is 
able to move in radial direction if instability indeed sets 
in. For Taylor vortices radial fluid motion is strongest 
near the gap center. This suggests choosing 



rp=ri + ^ A(a%). 



(6) 



Experimental data |18j indicate that this might be rea- 
sonable also for flow patterns other than Taylor vortices. 

Let us first evaluate the critical Reynolds number for 
resting outer cyhnder defined as Ri^c = Ri{R2 — 0)- For 
i?2 = we have Vn = r2, £ = ad and = (ri -I- ?'2)/2, 
yielding from (0) 



RiAv) = 



(l+7?)2 



2r,V(l-ry)(3 + r;)^ 
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where all constant factors meant by have been ab- 
sorbed in a. We fix a by matching Ri^dfj) to linear 
stability theory in the narrow gap limit, because in this 
limit Ri c is known very accurately, e.g. Ppp^. The 
combination 4(1 — r])Rl^/{l + rj) = is the critical 
Taylor number for i?2 = whose value for 77 ^ 1 is 
2/a'^ = TcKi 3416 §|. Thus we have a = 0.1556 corre- 
sponding to ^ = (i/6.427, a reasonable value, which was 
also found by the same reasoning in . 

Since Ri.c diverges for 77 — > 1 and 77 — > 0, one better 
rescales Ri and R2 to discuss the stability curve i?i(i?2) 
for general radius ratio 77, cf also We employ TZi = 
Ri/Ri.c a-ud TZ2 = R2/ {vRi.c)i which will happen to be 
useful also in the limit when the Reynolds numbers go to 
infinity. 

After eliminating a — £/d from (^) with ^ we arrive 
at the final formula for determining the stability bound- 
ary: 



(7^l-7^2)' 



n p 



=A(4)-^ (8) 



(l-77)(3 + 77) 



We now calculate a as explained above and get 



aiv) = (1 - 



I (l+7?)3 

2(1 + 377) 



(9) 



a increases weakly with 7; from w 1.4 to 1.6. We always 
took this 0(77), but differences to taking a — const within 
this range could hardly be noticed in our following fig- 
ures. Since ac?„ can also be understood as the radial 
size of a Taylor vortex for counterrotating cylinders, this 
range of a can be confirmed by inspecting several figures 
in Refs. |-||,|. 



1]-^ + 4t]+1 / 3 + 77 



2(1-1-77)2 V 1 + 377 



(10) 



is also weakly dependent on 77, ranging between —0.75 
and « —0.90. The continuous differentiability of the sta- 
bility curve 7?.i(7?.2) at adn/d — 1 can not trivially be 
satisfied. But it is possible because A enters twice in (H) , 
namely on the rhs. via £, cf (||) and (|^), and on the Ihs. 
via Tp, cf (H). The argument of A is the same in both in- 
stances. The curvature of the boundary line has opposite 
sign to the left and to the right of TZ'2 for all 77. 

Eq. d) together with (|), (|)-(0), and (|) constitute 
the main result of this work. It is an implicit equation 
determining the full stability curve TZi{TZ2) for all radius 
ratios 77. It is easily solvable; analytically in special cases, 
and numerically in general. 

Although we introduced several parameters, only a has 
been adjusted to linear stability theory. The others have 
been fixed by simple physical arguments and are not ar- 
bitrary fit parameters. 
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3. Comparison with data and stability theory 

We now discuss the contents of (||), i.e., how the sta- 
bility boundary globally looks like. 

Fig. |l| shows i?i,c(?7), Eq. (0), compared with Ri^dv) 
from the Navier-Stokes linear stability theory and some 
selected experimental data. The good overall agreement 
for rj approaching 1 and even for small 77 comes because 
Eq. (|^) has the correct asymptotic behavior, cf |^,|| for 
77 ^ 1, for 7/ 0, and The I/7; law may be 
considered as a consequence of our choice for rp, since 
with Vp = Ti [ p^ , for example, it is Ri^c = const for 77 — > 
0. Also the prefactor comes out remarkably well. We find 
V^i.c = 11-9 in the limit 7; — > 0, whereas rjRi^c ~ H from 
linear stability theory ||^, cf also the inset of Fig. |l|. Some 
previous authors gave phenomenological scaling formulae 
which, however, do not show the correct exponents of 
the asymptotic laws for 77 ^ or 7/ -> 1 @. 

In particular. Coles finds Ri^c oc l/{riy^\n l/rj) for 77 
which does not agree with linear stability theory 
(see the inset of Fig. 1 which clearly does not show a 
logarithmic singularity). The origin of this logarithmic 
singularity is that he took an integral mean of — L^(r)'/r^ 
over the range ri < r < r„ for R2 < 0, and ri < r < 
r2 for i?2 > 0, respectively, instead of choosing a well- 
localized position Tp as we do. 



In Figs. the corresponding comparison is offered 
for the stability boundary (|8|), in the original variables 
i?i(i?2), for 7? = 0.20, 0.50, and 0.964, covering the whole 
range for which experimental data seem to be available. 
Intermediate values of 77 as in Refs. P,pl]-jl3[| give essen- 
tially the same satisfying picture. 

For i?2 > we find excellent quantitative agreement, 
and for i?2 < there still is satisfactory agreement. The 
remaining discrepancies may be attributed to several rea- 
sons. First, the experimental data are affected by the fi- 
nite length of the cylinders which becomes important for 
sufhciently negative R2 j2^JT^. This can be clearly seen 
in Fig. ^ as a deviation of the Donnelly-Fultz data from 
linear stability theory, the latter being calculated with 
aspect ratio infinity. Snyder's data in Figs. p| and ^ are 
supposed to be free of finite length effects jfW . Still, the 
difference of his data to our curve in Fig. Hnear i?2 = 
is significantly larger than in other cases if we take the 
77 value reported in |l^], rj = 0.959. But the deviation 
is much below the claimed 1% error in the experimen- 
tal value for 77 iQ ; we indicated the corresponding error 
in i?i,c by a bar. Since the deviation could be caused by 
the strong 77 dependence ~ 1/ ^Jl — 77, cf Fig. |l], we match 
the i?i,c values of the experimental data and of the lin- 
ear stability theory (and thus our Eq. (^) by choosing 
rj — 0.964. Now we find the same satisfactory overall 
picture as for the other 77 values. 

Our curves correctly show a minimum at some nega- 
tive i?2. This can be attributed to the crossover in the 
boundary condition from solid to free, a = 1 leads to a 



non-differentiable corner-type minimum at i?2 = and 
to considerably higher curves for i?i(i?2) at counterrota- 
tion, as was also found in jl^]. Further on, taking a > 1 
and i?2 = as in fl^ gives a minimum at i?2 = which 
is not even continuous. 

It is not surprising that our minimum's position 



7^2 



1 



2(1 + 77) V 1 + 37? 



3 + 77 



its depth and width do not match the measurements too 
well since the minimum is located in the crossover region 



R*2 < R2 < 0. 



We have marked R2 with an arrow in 



the figures. Choosing some other, smoother crossover 
function A in (H) allows to modify the location and the 
width of the minimum. 

Interestingly, for 7; — > 1 the minimum is exactly at 
i?2 = 0, cf Fig. ^, independent of the value of a (as long 
as a > 1) and of the precise form of A, and nicelv com- 
patible with experiment j|jl^ and theory [§,||J|,^,^. 

For R2 far left of the minimum our curves calculated 
from (H) are clearly systematically too high. Before dis- 
cussing this, we complete our analysis by looking at in- 
teresting special cases, in which explicit formulae follow. 



4. Asymptotic formulae 

As long as the vortex flow fills the whole gap, i.e., for 
1^2 >'R'2^ we have A{adn/d) = 1; then Eq. (^ can be 
solved for TZi: 




3 + 77 

for 7^2 > n*2 



(11) 



This implies a leading correction cx 7?.^^ relative to the 
Rayleigh line TZi = TZ2 for large TZ2- 
For 7?.2 < 7^2 ) llie quantity 



^"Vi'' ^ 7/2 7^l-7^2 



(12) 



is a suitable parameter for discussing the strongly coun- 
terrotating case as well as the narrow and wide gap limits. 
For e — > 0, corresponding to r„ — > ri, we find from (H) 

-rj'n2 - (!^y^'-rj'nu 7^l » {l-rj^/'C^/^ (13) 



Civ) 



2^(3 + 77) 



1/5 



^a4(2-a)(l +77)7^ 
The leading term for large \Tl2\ reads 

7^l = c(77)|7727^2|'/^ 7^2^-oo. 



(14) 



The exponent 3/5 was already found by Donnelly and 
Fultz |iq] by a dimensional argument. That, of course. 
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did not yet give the complete 77 dependent prefactor 
rf/^C{r]) and also not the non-trivial corrections con- 
tained in (|T^). These turn out to be non-neglegible 
in the range of TZ2 accessible to experiment. Since 
£ oc \rf'TZ2\~'^^^ for sufficiently large |ry^7?.2|, the conver- 
gence to the leading behavior (13) and (|lj) is rather bad 
in terms of TZ2 for fixed 77 and becomes worse with de- 
creasing rj. 

Coles iQ has chosen the value of a by matching his 
analytic formula to linear stability theory for 77 ^ 1 in the 
limit i?2 ^ —00. This choice leads to an 77 independent 
value of a R:! 1.3. This means he adjusted his formula to 
two values from linear stability theory whereas we only 
use one. 

On the other hand, for 77 ^ 1 Eq. (|l3| ) becomes exact 



for aU 7^2 < 11% 
a = 8/5 from (10) and 



In this limit we have TZX 



-3/4 and 

Furthermore, the position of 
the minimum approaches 7^.2, min = 0, and the stabilit 
boundary becomes symmetric for \'R-2\ < |7^2l; cf Fig. ^ 
Therefore, in Fig. ^ the asymptotic formula (|l^ would 
be indistinguishable from the full curve, whereas in 
Figs. 1^ and ^ it would be a rather poor approximation. 

Of particular interest is the wide gap limit i] — > 0, since 
it has not yet been examined very much. It is accompa- 
nied by large e for any finite TZ2 interval, cf (|2|). For 
e ^ 00 the relevant expansion parameter in (H) is e^^^^. 
We keep the two leading terms and consider rj small in 
them: 



3 + 2^2^ 

1 + V2 ' 



V3 



7^2, 



(15) 



for - Ui/tj^ < 7^2 < 7^2 



Here 11^ ^ -\/3/2 and a = \/2 for 77 < 1, cf (|T^) and 
(^. Because of the bad convergence in terms of e"^^/^ to 
this limiting behavior the line Tii{TZ2) according to ( p^ ) 
is clearly visible only for rj < 0.01. Of course, for fixed 
77 > and sufficiently largely negative TZ2 equation ( p^ ) is 
valid again because then e ^ instead of — )■ 00. As can 
be seen from (14) and ( p^ ) the most reasonable variables 
for strongly counterrotating cylinders, i.e., TZ2 —00, 
are TZi and 7fTZ2- Note that the meaning of their ratio is 
77^7?,2/7?.i = u!2/ll)i. This frequency ratio as a parameter 
has already been used in 



Treating p as a free fit parameter we find the re- 
maining discrepancies between our stability boundary 
and that of experiment or linear stability calculations 
nearly removed for p = 0.47,0.46,0.43, corresponding to 
77 — 0.964,0.5,0.2, respectively. This is exemplified in 
Fig. ||by the dotted line; for the other 77 values the agree- 
ment is similarly good. Position shifts of this magnitude 
are confirmed by several figures in Refs. which 
show that the center of each Taylor vortex is situated 
nearer towards the inner cylinder. Analytically we find 
roughly a ^ 1/p and C approximately ~ p if p is near 
1/2 and for 77 not to small. We could not figure out a 
simple physical condition for fixing p(ri). For example, 
second order continuous differentiability of the stability 
boundary at TZi, turns out not to be realizable. 

Our argument does not take into consideration non- 
axisymmetric perturbations of the laminar flow. As we 
know from experiment and linear stability theory [^j , 
the relevant perturbations are indeed non-axisymmetric 
beyond a finite negative i?2- They lead to a somewhat 
lower stability boundary than that determined with ax- 
isymmetric disturbances. The difference for the values 
of i?2 considered here is pretty small. It turns out that 
our deviations from the correct stability boundary are of 
the same order of magnitude. Furthermore these devi- 
ations start already near i?2,min- We therefore think it 
to be acceptable that the possibility of non-axisymmetric 
disturbances is neglected. 

Of course, also other corrections have to be considered. 
For example, energy dissipation is curvature dependent. 
Therefore it seems natural to extend the rhs. of (^) to 
(i^^/^^)(5r)^(l -I- M/rp) with b a constant, i/vp vanishes 
for rj 1 or for TZ2 ^00 and is rather small even 
for 77 — > 0; for example in the case TZ2 = we have 
i/vp — a Pi 0.16. Thus ( p^ ) remains valid, but with 
a different a in C. Note, without this curvature cor- 
rection the value for 0(77), Eq. (||), is close to a = 8/5 
which minimizes C for fixed 77 but variable a. In general, 
this means the curvature correction shifts the stability 
boundary upwards if TZ2 is well below TZi, counteracting 
the influence of the rr, shift. 



6. Summary 



5. Improvements 

Now we discuss possible corrections to our main result. 
It turns out that the precise value of Vp has a surprisingly 
well observable influence on the stability curve. 

To see this, we generalize (0) to 



ri +pdA(a^) 
a 



(16) 



with the new parameter p describing the position of the 
disturbance; p = 1/2 gives the previous case, cf (|), to 
locate it in the center of the relevant gap. 



Similar arguments as are presented here for Taylor- 
Couette flow have been invoked to elucidate the phys- 
ical mechanism of the Rayleigh-Benard instability pi[ . 
But in the Taylor-Couette system it carries much fur- 
ther since this has two additional parameters, 77 and R2, 
whose influence on the shape of the stability boundary 
we could successfully explain. At Rayleigh-Benard in- 
stability a corresponds to the ratio of the a's for the 
two distinct cases of free and solid boundary conditions. 
From 1^ we flnd for this ratio w 1.13. We understand this 
value closer to 1 than we found here for the counterrotat- 
ing Taylor-Couette system from the stronger influence of 
the surface tension in an open Rayleigh-Benard system in 
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contrast to the completely missing surface tension at the 
nodal surface r„ in Taylor-Couette flow. Therefore this 
nodal surface is much softer, allowing for a larger exten- 
sion of the Taylor vortices beyond r„. One could check 
this by putting another fluid layer on top of the origi- 
nal Rayleigh-Benard fluid; if this diminishes the surface 
tension, a larger a should result. 

To conclude: While accurate results for the onset of the 
Taylor-Couette instability derived by linear stability the- 
ory are available, in this paper we have drawn a picture 
that gives comprehensive physical insight into the insta- 
bility mechanism and yields simple analytical expressions 
for the stability boundary in the full parameter space. 
These agree astonishingly well with experimental data 
and with the rigorous theory, especially in regard of how 
straightforward the argument is. The position where the 
initial perturbation occurs and the change of the type of 
boundary when the nodal surface moves inwards turned 
out to be important for the wide gap behavior and for 
the existence of a minimum of the stability boundary 
for negative R2- The values found for the parameters a 
and a are reasonable and agree with experiment as well. 
Assuming a radius ratio dependent small shift of the per- 
turbation's initial position towards the inner cylinder we 
were able to understand the remaining quantitative dis- 
crepancies. 

We believe that these simple ideas should be appli- 
cable to instabilities in more complicated rotating fluid 
flows, where the standard linear stability analysis is too 
difflcult. 
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FIG. 1. Critical Reynolds number Ri,c as a function of 
77 from Eq. (Q) (solid line) , from linear stability theory [^] 
(dashed line), and from various experiments |H,P_0l-h_3| (sym- 
bols). For the error bar see Fig. ^ Inset: The same for _Ri,c 
rescaled to Ri^c vV^ ~ V- 

FIG. 2. The stability line i?i (i?2i for rj = 0.2 from Eq. (|) 
(solid line), and from experiment [[13| (symbols). The dotted 
line is our result with the perturbation's position parameter 
set to p = 0.43, see Eq. (|lq). The arrow marks R2- Inset: 
The same on a magnified scale. 

FIG. 3. The stability line i?i {R2) for r7 = 0.5 from Eq. (|) 
(solid line), from linear stability theory [blM] (dashed line), 
and from experiment |l^ (symbols). The dashed line may be 
reproduced with p — 0.46. The arrow marks R2. Inset: The 
same on a magnified scale. 
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FIG. 4. The stability line i?i(7?2) for r) = 0.964 from 
Eq. (solid line), and from experiment (symbols). The 
vertical bar denotes the error in _Ri,c due to = 0.959 ± 1% 
[Q. p = 0.47 fits the experimental data. The arrow marks 
i?2- Inset: The same on a magnified scale. 



